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Abstract. Two derivations are given for a matrix- valued func- 
tion ^c{z), defined on the d-torus, that can be associated with 
a discrete, translationally periodic bar-joint framework C in M.'^. 
The rigid unit mode (RUM) spectrum of C is defined in terms of 
the phases of phase-periodic infinitesimal flexes and is identified in 
terms of the singular points of the function z — > rank {z) and 
also in terms of the wave vectors of excitations with vanishing en- 
ergy in the long wavelength limit. To a crystal framework C in 
Maxwell counting equilibrium we associate a unique multi- variable 
monic polynomial pc(zi, . . . , Zd) and for ideal zeolites the algebraic 
variety of zeros of pc{z) on the d-torus determines the RUM spec- 
trum. The matrix function is related to periodic "floppy modes" 
and their asymptotic order and an explicit formula is obtained for 
the number of periodic floppy modes for a given supercell. The 
crystal polynomial, RUM spectrum and the mode multiplicity are 
computed for a number of fundamental examples. In the case of 
certain zeolite frameworks in dimensions 2 and 3 direct proofs are 
give to show the maximal floppy mode property (order N). In par- 
ticular this is the case for the cubic symmetry sodalite framework 
and for some novel two-dimensional zeolite frameworks. 

1. Introduction 

Let C be a mathematical crystal framework, by which we mean a 
connected structure in the Euclidean space M*^ consisting of a set Cg of 
framework edges (bars or bonds) with a corresponding set of frame- 
work vertices (joints or atoms) such that Cg is periodic with respect to 
a discrete translation group T of isometries of M'^ with full rank. We 
mainly consider the case d = 3 together with the physical locally finite 
assumption in which Ce is generated by the translations of a finite set of 
edges and vertices. Such a geometric bar-joint framework C serves as a 
model for the essential geometry of the disposition of atoms and bonds 
in a material crystal Ai. In this case the vertices have atomic identi- 
fiers, such as H, He, Li, B, ... , and the chosen edges may correspond 
to the strong bonds. The identification of strongly bonded molecular 
units, such as Si04 and TiOg, imply a polyhedral net structure for C 
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and in particular aluminosilicate crystals and zeolites provide a fasci- 
nating diversity of tetrahedral nets in which every vertex is shared by 
two tetrahedra. 

Material scientists are interested in the manifestation and explana- 
tion of various forms of low energy motion and oscillation in materials. 
Of particular interest are the rigid unit modes (RUMs) in crystals and 
zeolites. These are the low energy (long wavelength) modes of crystal 
lattice oscillation which are related to the relative motions of rigid units, 
such as the Si04 tetrahedral units in quartz. The wave vectors of these 
vibrational modes are observed in neutron scattering experiments and 
have been shown to correlate closely with those for the modes observed 
in numerical simulations with periodic networks of rigid units. In both 
the experimental measurements and in the computer simulations the 
background mathematical model is classical lattice dynamics and the 
rigid unit modes are observed where phonon dispersion curves display 
markedly low energy. There is now a considerable body of literature 
tabulating the rigid unit modes of various crystals and it has become 
evident that the primary determinant for these low energy modes in 
a material A4 is the geometric structure of its abstract framework C. 
This work shows that typically the wave vectors of RUMs lie along lines 
and planes in reciprocal space but there are surprises. Firstly the wave 
vectors often lie on curved "exotic" surfaces in reciprocal space 
Secondly some zeolites have one or more rigid unit mode oscillation for 
every wave vector, representing an extraordinary degree of flexibility 
ig, H]. See also [5] and [21]. 

As we shall demonstrate the mathematics of rigid unit modes is 
essentially a linear first order theory and one can side-step much of the 
lattice dynamical formulations that relate to higher energy phonons 
and their dispersion curves. This was outlined in the seminal paper of 
Giddy, Dove, Pawley and Heine [3] and has been pursued effectively in 
the recent mathematical analysis of Wegner [22]. 

Motivated partly by material phenomena, in Owen and Power [16] 
the rigid unit mode spectrum was formally defined for an abstract 
crystal framework C with its translation group T as the set Q = Q{C) of 
points u = (wi, U2, W3) arising from the nonzero infinitesimal flexes u of 
C which are T-periodic modulo the phase factor oj = (wi, U2, 003). Here 
(x!j = exp(27rikj) and k = (fci, ^2, fcs), located in [0, 1) x [0, 1) x [0, 1), is 
the reduced wave vector associated with u. The rigid unit modes of C 
are defined to be the infinitesimal flexes u and so one may define the 
multiplicity of rigid unit modes for the pair (C, T) with wave vector k 
(or phase point u) as a vector space dimension. See Definitions 12.61 



The infinitesimal flex perspective is useful in that it brings into play 
the well-established theory of infinitesimal rigidity for finite bar-joint 
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frameworks. We show here that the consideration of general infinites- 
imal flexes in infinite bar-joint frameworks is effective in the analysis 
of the RUM spectrum. For example the identification of infinitesimal 
flexes localised to lines and planes (and appropriately partially phase- 
periodic) leads to planes and lines (respectively) in the RUM spectrum. 
We also see formally, in Theorem 15.11 that the vector space formula- 
tion of rigid unit modes for C does correspond to the phenomenological 
appearance of rigid unit modes in materials and in simulations. That 
is, the phase of a phase-periodic infinitesimal flex corresponds to a 
wave vector for which the harmonic oscillation excitations of vertices 
(atoms) have vanishing bond deformations in the infinite wavelength 
(slow oscillation) limit. 

In [in] various flexibility and rigidity properties were considered for 
abstract crystal frameworks in d dimensions. In the case of square- 
summable displacement vectors this lead to an associated matrix- valued 
function defined on the rf-torus T'^. This function, ^c{z), arises from 
a representation of the infinite rigidity matrix for C as a Hilbert space 
multiplication operator for a matrix- valued function, known in Toeplitz 
operator theory as the symbol function. The RUM spectrum Q{C) was 
shown to correspond to the set of points in T*^ where the null-space of 
$(^) is nontrivial. In the case of periodic tetrahedral nets (mathemat- 
ical zeolites) this is the set of zeros for the determinant of ^c{z) and 
the so the RUM spectrum of the crystal framework is identified with 
the zero set of a multi-variable trigonometric polynomial on T*^ with 
real coefficients. 

In what follows we give two direct derivations of $c(-2) which avoid 
Hilbert space rigidity operators and we formalise the association of a 
unique multivariable polynomial with a well-constrained crystal frame- 
work. The first derivation is elementary and arises from the identifi- 
cation of the family of matrices that characterise the rigid unit modes 
as null eigenvectors. The second arises from an explicit block diago- 
nalisation of the finite rigidity matrices associated with supercells and 
here ^c{z) appears as the "completion" of matrix- valued "functions" 
on finite subsets of the rf-torus. A companion viewpoint to this is that 
the supercell rigidity matrices associated with fixed period infinitesi- 
mal motions are obtained by "sampling" ^c{^) at appropriate rational 
points of the (i-torus. 

Using block diagonalisation we analyse the asymptotic order of the 
number of periodic infinitesimal flexes, or floppy modes, associated with 
supercells. Also we give an explicit formula for the number of periodic 
floppy modes associated with an ri x r2 x supercell, which resolves a 
question posed by Simon Guest. This number is computable from the 
matrix function as the sum 

dimker($c(i^'')) 

0<ki<ri,l<i<3 
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where k = {ki,k2,h) and = (gS^^^iAi, e^'^^^^Aa^ g27ra-3/r3)_ g^jg^ 
give a formal proof for a general crystal framework in M'^ that C is 
of "order A^" if and only if Q{C) has dimension d (and hence coincides 
withT*^). Here "order AT" means that with increasing supercell size the 
number of periodic floppy modes increases as the order of the number 
TV of atoms in the supercell. While it is generally true that in three 
dimensions, for example, order N"^/'^ (resp. order N^^^) is equivalent 
to RUM dimension 2 (resp. dimension 1) when Q{C) is standard, this 
may be in question if Q{C) is exotic and does not admit a rational 
parametrisation. We also prove, by straightforward linear algebra, that 
order is equivalent to the existence of localised floppy modes. 

In the final section we compute the rigid unit mode spectrum, the 
mode multiplicity and the crystal polynomial for various fundamen- 
tal examples including the frameworks for /3-cristobolite, sodalite and 
perovskite. We give the novelty of a two-dimensional framework with 
an exotic spectrum. This is defined by the simple requirement that 
the unit cell consist of an 8-ring of congruent equilateral triangles bor- 
dering a regular octagon. Its crystal polynomial pc{z,w,u) has total 
degree 6 and two exphcit cubic factors which contribute two "figure 
eight" curves to the spectrum. Also we indicate an approach to ex- 
otic curves and surfaces in terms of RUM spectrum evolutions in which 
the spectrum evolve in response to the "phase transition" through a 
periodicity-preserving framework deformation. 

With regard to order A^, we show that it is possible for a "two- 
dimensional zeolite" to have order A^, which resolves a problem posed 
by Mike Thorpe. In three dimensions it has been known from exper- 
imental studies that the cubic symmetry form of sodalite exhibits a 
rigid unit mode for each wave vector. See Dove-et-al [1], [6] where this 
mode of flexibility is linked to material properties. We provide a direct 
geometric proof that Csod has a localised infinitesimal flex and so for 
this reason has order A^. 

The focus in the sequel is for idealised crystal frameworks and the 
RUM spectrum and mode multiplicity comply with the strict mathe- 
matical sense of the definitions. In particular this spectrum is a real 
algebraic variety. On the other hand material crystals, such as the 
aluminosilicates surveyed in [1], have a rigid unit mode distribution 
determined experimentally. The accord between some of these sets of 
wave vectors underlines the efficacy of the rigid unit model as a primary 
determinant of low energy excitation modes in materials. 

The development here, which is self-contained, relates in part to 
themes and problems discussed at a London Mathematical Society 
workshop on The Rigidity of Frameworks and Applications. Q We 
briefly indicate some very recent papers that also address the rigid- 
ity and flexibility of periodic structures. On the mathematical side. 
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Borcea and Streinu pQ and Malestein and Theran jT3] examine peri- 
odic frameworks which are otherwise generic and in particular the latter 
paper obtains a combinatorial characterisation of periodic infinitesimal 
rigidity for frameworks in M^. Ross, Shulze and Whiteley [12], on the 
other hand, consider periodic frameworks which are rich in symmetry 
but are otherwise generic and obtain counting conditions for the pre- 
diction of finite motions. In Power [17] we derive the Borcea-Streinu 
rigidity matrix given in [1] and obtain Maxwell- Calladine and Fowler- 
Guest style equations (see [9], [15], [20]) for periodic frameworks. The 
flexibility considerations in these four papers are liberal in allowing 
affine deformation of the (Bravais) lattice structure while maintain- 
ing framework periodicity with respect to the deforming lattice. Such 
flexibility contrasts with the considerations here in which, on the one 
hand the lattice is fixed, an appropriate constraint for RUM analysis, 
and on the other hand periodicity for infinitesimal flexes is relaxed to 
periodicity-modulo-phase. Nevertheless there are interesting connec- 
tions to explore between these perspectives, such as the analysis of 
RUM evolutions indicated in Section 5. On the applications side, in 
addition to the work on RUMs indicated above we note that Kapko, 
Dawson, Treacy and Thorpe [12] have computed the flexibility win- 
dows (related to minimum and maximum volume deformation) of the 
frameworks of ideal zeolites. 

The author is grateful for stimulating discussions and communica- 
tions with Martin Dove, Simon Guest, John Owen, Mike Thorpe, Franz 
Wegner and Walter Whiteley. 

2. Crystal frameworks and infinitesimal rigidity. 

We first define crystal frameworks and their matrix-valued functions 
and give the main ingredients of infinitesimal flexibility for translation- 
ally periodic bar-joint frameworks. 

The following notational conventions are adopted. We write k for the 
label of a framework vertex in the unit cell, while k denotes a general 
multi- index (/ci, . . . , kd) in and k denotes a wave vector, as indicated 
above. 

2.1. Crystal frameworks and motifs. We focus mainly on three- 
dimensional frameworks since the adaptations needed for the extension 
to d = 2,4,5,... are essentially notational. Let G = {V,E) be a 
countable simple graph with V = {vi,V2, ■ ■ ■}, E C V x V and let 
Pi = p{vi) be corresponding points in M^. Then the pair {G,p) is said 
to be a bar-joint framework in with framework vertices, or joints. 
Pi and framework edges, or bars, given by the straightline segments 
[pi,Pj] between pi and pj when {vi,Vj) is an edge in E. 

An isometry of M.^ is a distance-preserving map T : — > M^. A full 
rank translation group T is a set of translation isometries {T^, : G Z^} 
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with Tk^i = Tk + Ti for all k,l, ^ I if k ^ and such that the three 
vectors 

a = Te,0, b = TejO, c = Te^O, 

associated with the generators ei = (1, 0, 0), 62 = (0, 1, 0), 63 = (0, 0, 1) 
of are not coplanar. 

Definition 2.1. A crystal framework C in M.^ with full rank translation 
group T and finite motif F^) is a countable bar-joint framework 
{G, p) in with a periodic labelling 

V = {v^^k : K e {I,. . .t},k e Z^}, p^^k = p{vK,k), 

such that 

(i) = {pk,o '■ 1^ ^ {1; • • • ^t}} and F^. <Z E is a. set of framework 
edges with at least one framework vertex in 

(ii) for each n and n the point p^.n is the translate T^p^.o? 

(iii) the set of all framework vertices is the disjoint union of the 
sets Tk{F^),k e Z3, 

(iv) the set Ce of all framework edges is the disjoint union of the sets 

r,(Fe),fcez3. 

The finiteness of the motif and the full rank of T ensure that Cy is 
a countable periodic set. We say that a crystal framework in is 
in Maxwell counting equilibrium if 3|-F„| = |-Fe|. Of particular interest 
amongst these are the 6-regular frameworks in which each atom is 
incident to 6 bonds and this includes the case of tetrahedral nets with 
each vertex shared by two tetrahedra. 

The principle ingredient of a crystal framework C is the pair [C^, C^) 
comprised by the periodic set Cg of framework edges together with its 
vertex set C^. Indeed one may view nonperiodic and quasicrystal frame- 
works in such terms. With translational periodicity there are choices 
for the isometry group T, since one may take subgroups, and there are 
further choices for motifs. It is natural to define the unit cell associated 
with T, and its labelling, as the closed parallelepiped set determined by 
the vectors a, 6, c. However it is also convenient, for notational book- 
keeping for the atoms on cell boundaries, to remove three of the faces 
of the unit cell to ensure no overlapping under translation. 

Definition 2.2. The (principal) partition unit cell for (C, T) is the set 
E of points with vectors tia + + ^sc, 0<tj<l,l<z< 3}, that is, 
the parallelepiped set determined by a, 6, c minus the three closed faces 
which do not contain a pair from {a, 6, c}. A general partition unit cell 
is an R^-translate of E. 

With E as in the definition it is natural to take F^ as the set C^P^E 
and to let F^. consist of the edges in E together with some framework 
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edges incident to framework vertices in and satisfying the partition- 
ing property (iv). Also it may be judicious, for symmetry considera- 
tions, to translate C to a "normalised" position so that the set Cy (1 E 
acquires certain simplicity or symmetry. 

Definition 2.3. A motif {Fy, F^,) is of incident type for C and S if 
Fy = Cy n E = {pk,o : I < k < t} and if each framework edge e in Fg 
has the form e = [pK,o,PT,s(e)] where k G Fy. The index (5(e) G is the 
exponent of e. 

Note that the exponent of e is (0, 0, 0) if and only if e, as a closed 
set, is contained in the unit cell E. Thus, the edges in Fe with nonzero 
exponent have one vertex in the partition unit cell and the other vertex 
in the T5(e) -translate of this cell. We may define the 

For example. Figure 1 shows such an incident type motif and parti- 
tion unit cell for the so-called kagome framework Ckag which is formed 
by corner-linked triangles in a regular hexagonal arrangement. The 
isometry group T in this case is generated by the two translations 
whose period vectors are given by the sides of the boundary of the unit 
cell. 



p 




Figure 1 . Motif and possible unit cell boundary for the 
kagome framework, Ckag- 



On the other hand Figure 2 illustrates a motif for the (rigid) squares 
grid Cgq made up of corner-linked rigid squares in regular square ar- 
rangement. This may be viewed as a 2d analogue of the basic octa- 
hedron network considered later which underlies the material crystal 
perovskite. The period vectors can be taken to be a = (2, 0),b = (0, 2), 
with Fy = {(0, 1), (1, 0)} for the motif vertices. Note that the edge 
64 = [p2,(o,i),Pi,{i,o)] in -^e = {ci, . . . ,65} has both vertices outside the 
partition unit cell [0, 2) x [0, 2). In particular this motif is not of incident 
type. 
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Figure 2. Unit cell and motif for the squares network, Csq. 

2.2. The function $c(^)- We now define the matrix- valued function, 
or symbol function, for C that is associated with a given motif [F^, Fg). 
Different motifs give matrix functions that are equivalent in the simple 
manner indicated below. 

Write z = {zi, . . . , z^), with Zi G C, |2;j| = 1,, to denote general 
points in the ci-torus. Typically u denotes a particular such point, 
especially a RUM value. Thus, with the usual multinomial convention 
UJ^ may denote the scalar value associated with a RUM point which 
is the product cu^' . . •cu^'*. On the other hand, just as z^ is convenient 
notation to denote the monomial function z ^ z^ on T, we write z^ 
for the monomial function z ^ z'' from T'^ to C. Since Zi~^ = "zl^ 
for points on the circle T we may think of general monomials z^ as 
products of the Zi or 'Zi with just non-negative powers. 

The examples we discuss occur in two and three dimension and in 
this case we dispense with subscripts and write {z,w) and {z,w,u) 
respectively for general points of and T'^ respectively. 

First we isolate the special case in which the motif is of incident 
type and has the added restriction that there are no "reflexive" edges 
of the form e = [pK,o,PK,<5(e)] in Fe- Note that the edge 65 in the motif of 
Figure 2 is of this form. Such simple incident type motifs are usually 
available for the frameworks of material crystals. 

The coordinates in of the point p^./c in will be written as 
{pK,i,k,PK,2,k,PK,3,k) and also, in the next section, as {pK,x,k,PK,y,k,PK,z,k)- 

Definition 2.4. Let C be a crystal framework with simple incident 
type motif (Ft,,Fe), with 

Fv = {Pk,o ■ I < i^<n}, Fe = {Ci-. 1 <i < m}. 

Then $c(^) is the m x 3n matrix of functions (pijiz) on T'^ such that 
for each edge ei = e = [pK,o,Pr,5{e)], 1 < i < m, 

0i,3(fc-l)+a(2;) = PK,afl " PT,a,5(e), 
(j)i^3l^r-l)+aiz) = -(Pk,<7,0 - PT,a,5{e))z^^''"\ 

with the remaining entries equal to zero. 
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The edges in a general motif have two forms, namely nonreflexive or 
reflexive. The former take the form e = [pK.,5{e,K),PT,s{e,T)], with k 7^ r, 
and with 5(e, k) and 5(e,r) possibly both nonzero. The "reflexive" 
edges have the form [pK.,5i{e),PK,52ie)]- We may choose such edges so 
that Si=0. 

Using a symbolic labelling for columns we now state the general form 
of ^c{z). 

Definition 2.5. Let C be a crystal framework in M.^ with general motif 
with 

Fv = {p^,o ■.l<K<n}, Fe = {ci: 1 <i < m}. 
Then $0(2) is the m x 3n matrix of functions such that for each edge 

e = [pK,5{e,K),PT,5{e,T)], with K T, 

{^C{z))e,i^,a) = {PK,a,5ie,K) - Pr,a,5(e,T))z^^''''^\ 

($c(2))e,(r,<7) = -{PK,a,5{e,K) - PT,a,S{e,T))z^^'''^\ 

while for each reflexive edge [pK,o,PK,5{e)] 

{^c{z))e,{K,a) = {PK,a,0-PK,a,S{e))i^ " 

with the remaining entries in each row equal to zero. 

Despite the notational detail above the algorithmic passage from 
crystal to matrix function is straightforward as we see in the following 
examples and in Section 5. Note that $(z) may be constructed directly 
from an associated standard motif graph G{Ai) which we define to 
be the vertex-labelled graph whose vertex set is V{Fy), the set of all 
vertices of the motif edges, whose edge set is F^, and whose labelling 
is given by asigning the k index to each external vertex. For example 
for the motif of Figure 2, the motif graph is also implied, with the two 
vertices of 64 carrying the labels (1, 0) and (0, 1). 

In the definition and notation above the dependence of $c(^) on the 
motif is straightforward. Replacement of a motif edge (resp. vertex) 
by a T-equivalent one results in multiplication of the appropriate row 
(resp. columns) by a monomial. Thus in general two motif matrix 
functions $(-2) and "^{z) satisfy 

^{z) = D^iz)A^z)BD2{z), 

where Di{z),D2{z) are diagonal matrix function with monomial func- 
tions on the diagonal and where A, B are permutation matrices, asso- 
ciated with edge and vertex relabelling. 

Example (a). Let Csq be the crystal framework in the plane defined 
by the motif illustrated in Figure 2. More precisely, with F„ equal 
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to the ordered set {(1/2, 0), (0, 1/2)} and = {ei,...,e5} 
connected network of rigid squares. It follows from Definition 12.51 that 



Csq is a 
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If the final row of $c^q(z, w) is deleted then one has the matrix function 
for the square grid framework, Cjp. say. This framework is in Maxwell 
counting equilibrium and 

det $Cg2 (^5 w) = A {w — z) {wz — 1) . 

Example (b) [16]. The matricial symbol function of the kagome frame- 
work Ckag determined by the six-edged motif of Figure 2 is given by 
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with determinant equal to a constant multiple of 
'zw{z — l){w — 1) {z — w). 

2.3. Polynomials for crystal frameworks. Let C be a crystal frame- 
work with a given isometry group T. If C is in Maxwell counting 
equilibrium, that is, if d\Fi,\ = \Fe\, then we may form the polyno- 
mial det($c(^)) of the matrix function associated with a particular 
motif. This is a polynomial in the coordinate functions Zi and their 
complex conjugates z^, that is, a multi- variable trigonometric polyno- 
mial on the d-torus. We remove dependence on the motif and formally 
define the crystal polynomial pc{z), associated with C and a lexico- 
graphic monomial ordering, as the product z'' det{^c{z)) where 7 is 
chosen so that 

(i) pc{z) is a linear combination of nonnegative power monomials. 



Pc[z) 



(ii) pc{z) has minimum total degree, and 

(iii) pc{z) has leading monomial with coefficient 1. 

Here the monomials are ordered lexicographically, so that, for exam- 
ple, the monomial function zfz2 has higher multi-degree than 2:1 ^f. It is 
thus possible to talk of the leading term of a multivariable polynomial 
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[3] . The crystal polynomial for the kagome framework, with motif and 
translation group as above is Pkag{z,w) = {z — 1) {w — 1) {z — w), or, 
expanded in monomial order, with precedence z > w > u, 

Pkag{z, W) = Z^W — Z^ — Zw'^ + Z + w'^ — W . 

We see in Section 5 that the factorisation of such polynomials is closely 
linked to the structure of the RUM spectrum of C. 

2.4. Rigidity matrices and floppy modes. Let C be a crystal frame- 
work in with framework vertices p^^k- A displacement vector for C is 
doubly-indexed sequence u of vectors UK,,k in regarded as a sequence 
of displacement directions, or instantaneous velocities, assigned to the 
frameworks vertices p^.fc- It is convenient to consider the vector space 
of all real displacement sequences, u = (uK^k), written 

Definition 2.6. A real infinitesimal flex u for C is a displacement 
vector u in T-iatom such that 

(^i^.fe - Ur,l,PK,k-Pr,l) = 

for each framework edge [pK,k,PT,i]- 

Thus a infinitesimal flex of C is an instantaneous atom velocity se- 
quence such that at the two ends of each edge the velocity components 
in the edge direction are equal. The terminology here, borrowed from 
finite framework theory, reflects the equivalent formulation which re- 
quires that, for all small time t the perturbed edge 

[PK,k + tU^^k,Pr,l + tUr^l] 

has length of second order smallness. That is, 

\[PK,k + tU^,^k,PT,l + tUr,l\\ = \pK,k-PT,l\+0{t'^) 

as t tends to zero. Here the implicit constant for O(t^) depends only on 
the edge and so, for example, an infinitesimal rotation qualifies as an 
infinitesimal flex (despite the unbounded nature of the flex sequence.) 

The rigidity matrix R{C) is an infinite matrix defined exactly as in 
the finite framework case. 

Definition 2.7. The rigidity matrix R{C) of the crystal framework C 
has rows labelled by the edges e = [pK.,k,PT,i], columns labelled by the 
framework vertex coordinate variables {k, x, k), [k, y, k), [k, z, k), and 
row e takes the form 

[■ ■ ■ (pK,fc - Pr,l) ■ ■ ■ {p-r,l - Pn,k) " " " ] 

where the entry {j)n,k —Pt,i) indicates that the three coordinates of this 
vector lie in the columns for (k, x, k), (k, y, k), (k, z, k). 
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We remark that one may take the view that R{C) is 1/2 J(C) where 
J(C) is the infinite Jacobian, evaluated at the p^^^k for the quadratic 
equation system 

\qK,k-qT,i\ =4 

where de is the length of the edge e. 

It is natural to take a linear algebraic perspective and consider var- 
ious linear transformations that derive from R{C). 

Let 

be the space of real sequences w = {we^k)eeFe,kez^ labelled by the frame- 
work edges. Then R{C) gives a linear transformation R : T-Latom ~^ 
Hbond- Indeed, each row of R has at most 6 nonzero entries and the 
image R{u) is given by the well-defined matrix multiplication R{C)u. 
As for finite frameworks one has the following elementary proposition. 

Proposition 2.8. The infinitesimal flexes of C are the displacement 
vectors which lie in the nullspace of linear transformation R. 

For r = (ri, r2, r^) in 1? let "H^ be the finite-dimensional subspace of 
displacement vectors u for which UK,n+m = u^^n if divides rrii,! < i < 
3. These vectors can be viewed as the displacements determined by 
their restriction to the r-fold motif, denoted {F^, FJ), for the indexing 
subgroup rZ^ = riZ x x r^Z. It is natural to take (resp. FJ) 
to be the union of the disjoint /c-indexed translates of Fy (resp. Fg), 
where k = (fci, /c2; ^s); < /cj < Tj — 1. Similarly define HI and note 
that R{C) provides a linear transformation 

of finite-dimensional vector spaces. For example this follows from the 
commutation relations in Section 3. 

In particular, for r = (1, 1, 1) the transformation R^"^^ is a linear 
transformation between the spaces of strictly periodic vectors. With 
respect to the natural basis this transformation provides a rigidity ma- 
trix which we denote as R{Ai) and refer to as the motif rigidity matrix 

m, m- 

The following definition accords with the usual usage in applications, 
where a crystal framework is said to be of order if the number of 
periodic fioppy modes possessed by a supercell is of the same order as 
the total degrees of freedom of the vertex atoms in the supercell. 

Definition 2.9. Let C be a crystal framework in and let R^'^^ be 
the finite-dimensional linear transformation induced by R{C) for the 
n-fold supercell, where r = {n,n,n). 

(i) A floppy mode or, more precisely, a periodic fioppy mode, of C is 
a nonzero vector u in the nullspace of R^^^ for some r G Z^. 
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(ii) A local floppy mode is a nonzero vector u = (ui^^k) in the nullspace 
of R{C) which vanishes outside some supercell, that is, Mk,^ = for all 
but finitely many values k. 

(ii) C is said to be of order N"" for floppy modes, where a = 0, |, | 
or 1, if for some constant C > the nullity z/„ = dimkeri?^^) satisfies 
Un > Cn^^ for all n, while there is no such constant for the power a + 1. 
In particular C is said to be of order if > cn^ for some constant 
c> 0. 

Note that if C has a local floppy mode u then by considering dis- 
jointly supported translates of u one may obtain linearly independent 
local floppy modes of order = ra^ in number in an ra-fold supercell. 
Further, these local floppy modes have n-fold periodic floppy modes 
associated with them and so it follows that C is of order A^ in this case. 

We now show that the converse implication follows from simple linear 
algebra 

Let us specify notation for the natural "basic sequences" of 'Hatom 
and Tibond- Write C,x,^y,^z for the standard coordinate basis of M^, 
C,x = (1,0,0) etc., and for a G {x,y,z} write iK,,a,k for the position 
indicator vector u in 'Hatom with 

where (5k,k' is the Kronecker delta. While T-Latom does not have countable 
vector space dimension its subspace of finitely nonzero sequences has 
the set {^K.<7,k} as a vector space basis. Similarly we write rj^^k for 
the basic sequence in l-Lhond which is zero but for the value 1 for the 
coordinate position e, k. 

For the subspaces of r-fold periodic vectors we introduce some anal- 
ogous notation. Write A; G r to denote an integer points k = {ki, /c2, k^) 
in with < /cj < rj — 1. Consider the set 

{L,a,k ■■ {1,2, . . . , |F^|},o- G {x,y,z},k er}, 

where ^K,a,k is the r-fold periodic vector whose k,', a', coordinate is 
SK,K'Sa,a'- This is a vector space basis for "H^tom- Also we write 

{r]e,k : e G {1, 2, . . . , |Fe|}, A; G r}. 

for the natural basis of T-ilond- 

Define the boundary Fj'*'^ of the motif (F^, Fg) as the set of vertices 
V{F(,)\F^. For r = {n,n,n) the boundary of the motif {F^,F^) is 
contained in the union of the boundaries of the motifs for the cells 
in the boundary of the n-fold supercell. It consists of the points that 
lie outside the union of the partition unit cells and in particular 
|^^r-j6drj/| ^ somc coustaut C and all n. In the simple case 

of a rectangular unit cell and an edge motif Fe of adjacent connection 
type, in the sense that the exponent 6^ for e in F^, takes one of the 
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seven forms (0, 0, 0), (±1, 0, 0), (0, ±1, 0), (0, 0, ±1), one can verify that 
|(^^r^Mry| ^ 3n'^\F^\ for all r = (n,n,n). 

Proposition 2.10. Let C be a crystal framework in M.^ which is of 
order N for floppy modes. Then C has a local floppy mode. 

Proof. With r = {n,n,n) let P be the linear "boundary projection" 
transformation on such that P^K,a,k = in,a,k if the pair k, k label a 
boundary point of FJ in and is zero otherwise. The dimension of 
the "boundary vector space" P{Wg) is equal to KFJ)^*^^^! and so is of 
order n^. 

Let M^, . . . , u^, with k = Un, be a basis for the nuUspace of R{Cn) = 
R^^\ Using the assumption, choose n so that k > dimP('HJ^). Thus 
the k vectors Pu^, . . . , Pu'^ are linearly dependent and for some set of 
coefficients Aj, not all of which are zero, we have 

= AiPn^ + ■ ■ ■ + AfcPn^ = P{Xiu^ + ■■■ + XkU^). 

Also, since u^, . . . ,u'' is a basis the vector u = Xiu^ + ■ ■ ■ + Afcu'^ is 
a nonzero vector in the periodic vector space with Pu = 0. Note 
that such a non zero periodic displacement vector, which is zero on 
the boundary vertices of the supercell, can be localised by truncation 
to create the floppy mode u', where m'^ q = ""k.o and u'^/^ = if k 
(0, 0, 0). To see explicitly that this truncation is an infinitesimal fiex one 
need only consider the (new) displacements, u'^.u'^ for the two vertices 
Pa,Pi3 of an edge e which has just one vertex, say, in the partition 
supercell. But in this case the internal vertex has a T^-translate to a 
boundary vertex of and so by periodicity = U/j = 0. Thus with 
the truncation change — = one still has the infinitesimal fiex 
requirement 

« - Up,pa-Pp) = 

□ 

3. The role of ^ciz). 

We first show how $c(-2) arises as a family of matrices parametrised 
by points z in the d-torus where the matrix for z = determines the 
infinitesimal fiexes that are periodic modulo the phase ZJ. Once again, 
for notational simplicity we assume that d = 3. 

3.1. Phase-periodic infinitesimal fiexes. Let ICatomi^hond be the 
complex scalar vector spaces ICatom = C (g) Hatom, JCbond = C (g) Hbond- 
Write /C^ for the complex vector space of displacement vectors v = 
(^K.fc) such that WK,fc = <^^Vk,o for k G Py, G Z^. This is a finite- 
dimensional subspace of JCatom and we similarly define /C^ C JCbond- 

The matrix R{C) provides a linear transformationn : /C^ — /C^ 
by restriction. Indeed, let 7,, 1 < i < 3, denote the usual generators for 
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1? and let Wi and Ui be the shift transformations on Katam and /Cbond 
respectively, with 

Then we have the commutation relations 

R{C)Wi = UiRiC), 1 < z < 3, 

and the identities WiU = uJiU, for u G /C^, and f/^t; = ZJit", for v G /C^. 
Thus for u in /C^, 

f/i(i?(C)M) = R{C){Wiu) = R{C){uJlu) = ulR{C)u 

and so R{C)u G /C^ 

Let {^re,cr : K e -Fi;, 0" G {x, y, z}} be the natural basis for the column 
vector space C^'''^"'. Write for the displacement vectors in /CJJ' which 
"extend" the basis elements ^k.o-.Oi ^ F^^o' G {x, Formally, in 

terms of Kronecker delta symbol, we have 

Similarly let r/e? e G -Fe, be the standard basis for C''^'=' and write 77^, e G 
Fg, for the natural associated basis for /C^, with 

Proposition 3.1. Let C he a crystal framework with matrix function 
^c{z) and letu G T*^. Then the scalar matrix $c(^) is the representing 
matrix for the linear transformation R^ : IQ^ — t- /C^ with respect to the 
natural bases {^^o-} ^'^^ {vt}- 

Proof. Let ( , ) denote the scalar product in /C^ for which the natu- 
ral basis elements are orthonormal, that is, using the Kronecker delta 
symbol, (fJeiVf) = ^ej- Then the entry of the representing matrix for 
R'^ in the e row and the k', a column is 

iRne,(.',.) = {RiC)C,.,Ve). 

From the definition of the rigidity matrix it follows that R^'^, ^ is a. 
sum of phase-periodic vectors provided by the edges e in Fg, with e = 
[PK,Si^Pr,S2]y where k' is equal to one, or both, of k and r. In fact, 
using the Kronecker symbol and writing Ve for the edge vector 

PkM -Pr,S2, we have 

To see this note that 

Also R{C)^'^i is cj-periodic and the vectors ?7g form a basis for the 

space of cu-periodic vectors in the range. Thus R{C)^'^, ^, and hence 
R^^k' ct' is equal to the sum indicated above. 
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Noting the form of $c(i^) in Definition 12.51 it follows 

□ 

The proposition gives a method for finding the RUM spectrum in 
the sense of the following preliminary definition. 

Definition 3.2. Let C be a crystal framework in M.^ with translation 
group T. Then the RUM spectrum Q{C) for (C, T) is the set of points 
u in T'^ for which there exists a nonzero vector in the nullspace of R^. 

The formalism here is essentially the one adopted by Wegner |22] 
in his computational determination of exotic RUM spectra. However 
here we bring the linear space transformations and structures more 
explicitly into play. We note also that Hutchinson and Fleck [11] in 
their engineering analysis of kagome related structures adopt a phase- 
periodic flex and stresses analysis, drawing on the formalism of lattice 
dynamics and Bloch's theorem. They also make explicit the kagome 
function matrix. 

3.2. $c(-2) and floppy mode counting. We now start afresh and 
derive $c(^) as the completion of a function defined on the rational 
points u = {wi, W2, w^) with Wi = e^'^*"' and a a rational number in the 
interval [0,1). Despite the notational density the proof is essentially 
the matricial form of a standard block diagonalisation argument for 
transformations that intertwine cyclic shifts. 

Let r = (ri,r2,r3) and consider the finite-dimensional space /C^ = 
7/J^®C in JCatom = Hatom®'^ cousistiug of r-periodic complex displace- 
ment vectors. This has a basis of vectors {^K,cr,fc}; as given in Section 2 
for T-Ua which we view as an orthonormal basis for a Euclidean (Hilbert 
space) norm. Similarly, /C^ is the space of r-fold periodic vectors in 
/Cfoond with orthonormal basis {?7e,fc}eeFe,fcer- 

Let oji = e^'^*/^', for / = 1, 2, 3, let be the finite discrete torus 

and consider the vector space of function on T'f. taking values in the 
column vector space C'^'^"'. This space has the same dimension as 
/C^ since, for example, one could take a basis consisting of the basic 
functions 

w -> Xfc(w)^K,<T, hie F^,a e {x, y, z}, k e r, 

where Xfc • C is the indicator function for the point u'' in T^. 

Here the set of vectors {^k,o-} is the basis for C'^''^"'. It is convenient to 
introduce a natural Hilbert space inner product so that the indicator 
functions above are orthogonal with norms {rir2r-j)~^^'^. This corre- 
sponds to viewing the function space as L^(T^) C'^'-^"' where L^(T^) 
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has Hilbert space (Euclidean) norm determined by normalised counting 
measure on T^. 

Similarly we introduce the function Hilbert space L^(T^) ® C'^=' 
which has the same vector space dimension as JCl- 

Note now that the vectors w'^ ® ^K,a, for k & r, given explicitly as 
the functions 

form an orthonormal basis. This is because of the standard orthogo- 
nality of the set of vectors 

(l,/3'',/32'',...,/3(P-i)'?), l<g<p, 

where (3 is the p^^ root of unity. 

Thus we may define the unitary transformation 

:/c:;^l2(t3)®c3|^"| 

in terms of the bijection 

^K,<7,fc w'' <S) ^K,<x, K e F^,,a e {x,y,z},k e r. 

In a similar way, for the edge motif vector space there is a natural 
Hilbert space isomorphism 

determined by a bijection between bases, 

?7e,fc w'' ® ?7e, e e Fe, k & T. 

Let \E'(w) be any complex function on whose values are complex 
matrices of size |Fe| x 3|F„|. Then this function determines a linear 
tr ansf ormat ion 

by pointwise matrix multiplication; for a function F{w) in the domain 
space {M^F){w) = '^{w)F{w). Note that one may view the transfor- 
mation as being in block diagonal form since it is evidently equivalent 
to a block diagonal transformation with matrix 

Lu eTf. 

The next proposition shows the role of $c(^) and provides a conve- 
nient representation of the rigidity matrices for supercell periodicity. 
In particular it provides a counting formula for the number of floppy 
modes associated with a supercell. 

Theorem 3.3. Let r G Zj^ and let ^^{w) be the function on the finite 
3-torus determined by the restriction of^ci^)- Then 

(i) the rigidity matrix transformation R^'^^ : /C^ — j- /C^ for r -periodic 
displacements is unitarily equivalent to the multiplication transforma- 
tion M^r, specifically, R^''^ = U^Mi^rUa- 
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(ii) the number of r -periodic floppy modes is 

dim ker i?'^''^ = dimker A^<j,r = ^ dimker($c(c<j^')) 

ker 

where u = (exp(27r/ri, . . . , exp(27r/ri) and the sum extends over the 
^i^2^z points of the finite torus.. 

Proof. Let k' ^ Fy,a & {x, y, z} and let k' E r and consider the expan- 
sion 

R'^'^^in' ,u,k' = -R(C)^K',o-,fe' = ^ <^(e,k),(K',a,k')Ve,k 

eeFe,k& 

in terms of the basis for /C|J. Thus the coefficients of the vectors ?7e,fc 
give the representing matrix for R^^\ with rows labelled by the {e,k) 
and columns labelled by the triples (k',(T, /c'). Let e = [Pk,Si,Pt,52] be 
an nonreflexive edge in Fg with edge vector Vg = Pk,5i — Pt,S2 ■ Then 
the translated edges TfcC, for k E r, have the same edge vector and the 
entries of the rigidity matrix for row (e, k) imply that 

tt(e,fc),(K,cr,fc+<5i) = (^'e)o-, 
«(e,fc),(T,<T,fc+<52) = —{Ve)a- 

Here k + 6i and k + 6i are taken modulo (the 3-tuple) r. All other 
coefficients for row (e, k) are equal to zero. 

With respect to the natural inner product on /C^ we have 

{R^'''^L,a,k+Si,Ve,k) = {Ve)a, 
{R^'^^iT,cr,k+52-,Ve,k) = —{Ve)a, 

with other entries of the row equal to zero. Thus, in view of the speci- 
fication of Ua and Ub we have 

{UhR^''>U:{w''+'' ® e.,.), W'' ® Ve) = -{Ve)a. 

With T = UbR^^'^U* this gives 

The case of a reflexive edge is the same except that there is at most 
one nonzero entry in the corresponding row, taking the form of the sum 
of the two values, namely {oJ^^ —uj^'^){ve)a. 

One the other hand, the form of $c given in Definition 12.51 shows for 
nonreflexive edges 

with a corresponding formula for a reflexive edge. Thus i?'-''^ and 
U^MijirUa have the same representing matrices and so are equal and 
the proof of (i) and (ii) is complete. □ 
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3.3. $c(-2) as a Hilbert space operator. Write /C^ and /C^ for the 
Hilbert spaces of square summable sequences in JCatom and JCtond- Then 
R{C) determines a bounded Hilbert space operator from /C^ to which 
intertwines the shift operators as before, although now these operators 
are unitary operators. Identifying square- summable sequences with 
square integrable functions one obtains unitary equivalences between 
K-l and L'^{T'^) ® C^l^-I and between ICf and L'^{T'^) ® C\^^\. With 
careful spatial identifications one can identify the corresponding (uni- 
tary) transform of the linear transformation R{C) with a multiplication 
operator between these matrix-valued function spaces where the mul- 
tiplying function is $c(^)- For more details see Owen and Power 
where other operator-theoretic considerations are indicated. 

We remark that the choice of z rather than z in the Definitions 
12. 4[ 12.51 for ^{z) is somewhat arbitary from the perspective of trans- 
lationally periodic crystal frameworks. However elsewhere we consider 
semi-infinite crystal frameworks and in this case it is natural to model 
displacement velocity sequences in the domain by analytic functions 
and power series, with positive powers of the variable. With this choice 
the negative exponent formalism of these definitions is implied and so 
it seems appropriate to build this in to the definition of ^{z) at the 
outset. 

4. Rigid unit modes. 

Let C = {Fy, Fe, T) be a crystal framework in M*^, = 2, 3, . . . , with 
motif {Fvi Fe), full rank translation group T and rigidity matrix R{C). 

Definition 4.1. (i) The RUM spectrum of C is the set Q{C) of points 
uj = [ui, . . . ,uJd) in T'^ such that there is a nonzero vector u with 
R{C)u = which is periodic modulo the phase u; 

0(C) = {w G T'^ : keri?^ ^ {0}}. 

(ii) The mode multiplicity function defined on Q{C) is given by 
/i(aj) = dimkeri?"^. 

(iii) The RUM dimension dimjiuMiC) of C is the real dimension of 
the real algebraic variety Q{C). 

It follows that 

n{C) = e : rank$c(u;) < d\F^\} 

and in the case of frameworks in Maxwell counting equilibrium, 

n{C) = {w G T"^ : det $c(w) = 0}. 

From the definitions we see that the periodic floppy modes, as we have 
defined them, are vectors in the nuUspace of R'^ for the points u in the 
rational ci-torus Tf^^ C T*^. Also we may consider the corresponding 
inclusion 
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From the remarks at the end of Section 2.1 it follows that the RUM 
spectrum and the rational subset fl/iop may depend on the choice of 
translation group, or, equivalently, on the choice of unit cell although in 
the notation we usually suppress this. However, we note the following 
simple relationship between the RUM spectrum Q{C,S) and the spec- 
trum Q{C, T) when 5 is a subgroup of T corresponding to a supercell 
of ri X ■ ■ ■ X Trf unit cells. If 

S = {Tfc : Tj divides ki,l < i < d} 

then there is a natural covering map 7 : T'^ — T'^ given by '~f{w) = 

' "' . It follows from the definition of RUMs that 

fi(C,5)=7(f^(C,r)). 

In particular the dimension of VL{C) is a property of C alone since the 
covering map is finite-to-one. 

The following theorem and the ensuing discussion relate the order of 
floppy modes to the RUM dimension of C. 

Theorem 4.2. (i) Let {C,T) be a d- dimensional crystal framework 
with RUM spectrum Q{C) C T'^. Then the nullity Un = dimkeri?„ of 
the rigidity matrix for n-fold periodicity satisfies 

d-i + \n{C) n T^l < z/„ < d\F^\\n{C) n 

where \F^\ is the number of vertices in the partition unit cell and where 
is the discrete torus for r = {n,n,n). 

(a) If Un > en"" for some c > 0,a > 0, then dim(n(C)) > a. 
Proof, (i) From Theorem 13.31 above we have 

Un = dimker/2„ = dimker($c('^^))- 

fcSr 

Thus the second inequality follows readily since dimker($c(cj) < d\Fy\ 
for all CO. On the other hand, if to'' E f2(C)nT^ then dimker > 
while for wave vector k = (0, 0, 0) we have dimker($c(l) •••!!)) ^ d, 
since there are certainly d linearly independent translation infinitesimal 
flexes. Thus the first inequality follows. 

(ii) follows from (i) since for any algebraic variety Q, if the dimension 
is less than the integer a then the cardinality of ^2 fl is at most of 
order □ 

Three illustrations will be given of the following corollary in Section 
5, including the Roman tile framework Cptom and the sodalite framework 

Corollary 4.3. With the notation above the following assertions are 
equivalent for a crystal framework C in W^. 

(i) C has a local floppy mode. 

(ii) > cn*^, for all n, for some c > 0. 
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(ill) dim(l](C)) = d. 
(iv) n{C) = T'^. 

Proof. The equivalence of (i) and (ii) follows from Proposition 12.101 
The last theorem shows that (ii) and (iii) are equivalent, and (iii) is 
equivalent to (iv) since Q{C) is a real algebraic variety in T*^. □ 

Remark 4.4. Suppose that Q{C) has dimension t with t less that d 
and that fl{C) is rationally parametrised in the sense that there are 
rational polynomials gi{s), . . . , gd{s), in variable Si, . . . , Sf , such that 

n{C) = {(exp(27rz^7i(s)), . . . , exp{2ingd{s))) : s E [0, l^} C T"^. 

Then it follows that Qfiop{C) is dense in Q{C) and that the n-fold pe- 
riodic floppy modes, which are associated with the rational points, are 
of order n*. In general however, in the case of exotic RUM spectra, one 
should expect examples where the rational points of the RUM spectrum 
are not dense. 

5. Crystal Frameworks. 

We first show how the points of the RUM spectrum Q{C) arise as 
the wave vectors k of the wave excitations of C which induce vanishing 
bond distortion in their long wavelength limits. This is followed by 
the determination of the RUM spectra, both standard and exotic, for a 
number of key frameworks. The infinitesimal flex approach here gives 
theoretical insights as well as alternative computational avenues and 
we intend to develop this further elsewhere. 

5.1. RUMs and long wavelength limits. Suppose that C is a crys- 
tal framework in M"^, with motif data {Fy, F^, T) and suppose that the 
vertices pK,k undergo a standard wave motion, 

VK,k{^) = PK,k + u^^kit), K E Fy,k eZ"^, 

where u^^kit) represents the local oscillatory motion of atom k in the 
translated unit cell with label k E Ij^. Following formula-simplifying 
convention, these motions are complex-valued, the case of real motion 
being recoverable from real an imaginary parts. (See Dove [8].) Thus 
we have 

UK,k{t) = MKexp(27rzk ■ k) exp(zat) 
where u = {Uf^)f^^F^ is a fixed vector in C"^'^"', where k is the wave 
vector and where a is the frequency. 

Considering the distortion A^it) for the edge e = [pK,,k,PT,k+5{e)] we 
have 

Ae(t) := bK,fc(t) - Pr,k+5{e)it)\'^ " ^^.^(O) - p^,fc+5(e) (O) ^ 
= 2i?e(pK,fc - PT,k+5{e),Uf,^k{t) - Ur,k+5{e){t)) 
+2Re{p^^k - PT,k+5(e),UK,k{^) - Mr,fc+<5(e)(0)) 

-|-e(u, k, /c, at) 
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where 

e{u,k,k,at) = \u^^k{t) - Uk+s{e){t)\'^ - kK,fc(0) - nfc+5(e)(0)p. 

First note that in any finite time period [0,T] the difference quantities 
e(u, k, /c, at) tends to zero uniformly, for all t G [0,T] and all k in Z^, 
as the frequency a tends to zero. This follows readily from the fact 
that for any 6 the quantity 

I sm{at + 6)- sin(at)p - | sin(^) - sin(0)p 

tends to zero uniformly for t G [0,T] as a tends to zero. 
For the other terms for Ae(t) note that 

which is zero, irrespective of t, if {u^u,^) is an infinitesimal fiex of the 
framework. 

It follows that we have proven the implication (i) implies (ii) in the 
following proposition. The converse assertion follows from the same 
equations and so we conclude that the points in f2(C) correspond to 
the wave vectors of RUM phonons that appear in simulations. 

Theorem 5.1. LetC he a crystal framework, with specified periodicity, 
and let k be a wave vector with point w e T^. Then the following 
assertions are equivalent. 

(i) {u''uK)K,k is a nonzero periodic-modulo-phase infinitesimal (com- 
plex) flex for C. 

(ii) For the vertex wave motion 

PK,k(t) = PK,k + u^exp(27rz k ■ k) exp(zat), 
and a given time interval, t G [0,T], the bond length changes 

^e(^) = IP^A't) - PT,k+5(e){t) \ - IPkA^) -Pr,fe+5(e)(0)|, 

tend to zero uniformly, in t and e, as the wavelength 271 /a tends to 
infinity. 

In the last two decades the RUM spectra of frameworks associ- 
ated with specific material crystals have been computed using the pro- 
gramme CRUSH [2]. The method here refiects principle (ii) in the 
proposition above and RUM identifications are obtained by lattice dy- 
namics simulation. This is based upon a double limiting process known 
as the split atom method. Here each shared vertex (often an oxygen 
atom) is duplicated, for each rigid unit, and connected by bonds of zero 
length and increasing strength, tending to infinity. In this setting the 
RUM wave vectors coincide with those for which the long wavelength 
limits have vanishing energy and so can be identified and counted. 
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Some of the results of this approach can be found in [3], [B], [T], 



5.2. Standard and exotic RUM spectra. The RUM spectrum of a 
simple crystal framework (C, T) in three dimensions often contains or 
consists of a union of a finite number of points, lines and planes of wave 
vectors k. In this case and in the order case, with Q{C) = T^, we 
say that Q{C) is standard. Otherwise we say that the RUM spectrum 
is exotic. (In view of the arithmetic subtlety alluded to in Remark 14.41 
one should view these terms as working definitions adequate for our 
discussion.) We shall see, for example, that the following basic crystal 
frameworks have standard RUM spectrum: the grid frameworks Cid 
in W^, the kagome framework, the kagome net framework in M^, and 
the octahedron net framework Cqct associated with perovskites and 
which may be viewed as a three-dimensional analogue of the squares 
framework Cgq. 

The infinitesimal flex perspective that we have been pursuing can 
enable the prediction of lines and planes in the RUM spectrum, even 
without the computation of the rank degeneracies of ^c{z), as we now 
indicate. 

Consider, for example, the basic grid framework C^a in the plane, 
as indicated in Section 2 but re-orientated so that the vertices are the 
points with integer coordinates and the translation group is that for 
integer translations. (The motif has a single vertex, = {p^}, and 
two edges.) Examining all edges it becomes evident that there exists 
an infinitesimal flex u supported on the x-axis, with MK,{fci,o) = (IjO) 
for all ki G Z. Using all the parallel translates Ti^q^i^^-^u of u, we may 
define the periodic-modulo-phase displacement v given by 



where 002 is a fixed point in T. Then v is well-defined and is a nonzero 
infinitesimal flex and so its phase {l,u}2) in lies in the RUM spec- 
trum. In the language of wave vectors we see that the RUM spectrum 
contains the line of wave vectors (0,a). By symmetry the line (a, 0) 
is also included. Similar arguments apply to the kagome lattice which 
also has linearly localised infinitesimal flexes. 

There is a general principle at work here which holds, as we illustrate 
below, for a crystal framework C which has an infinitesimal flex u which 
is 

(a) band limited, in the sense of being supported by a set of framework 
vertices within a finite distance of a direction axis for T, and 

(b) periodic for the direction axis direction. 
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Indeed, by (a) one can form a sum analogous to that above, using 
the complementary direction axis, to obtain a well-defined periodic- 
modulo-phase infinitesimal displacement. By linearity, and commuta- 
tion of the rigidity matrix with translations, v is an infinitesimal fiex. 
Similarly, if (b) is changed to the requirement of unidirectional phase- 
periodicity for phase oji then we deduce that cui x T is contained in the 
RUM spectrum. 

Similarly one can construct phase-periodic fiexes for the grid frame- 
work in three dimensions and higher. For the grid framework Cj2 in 
three dimensions one deduces from the evident line-localised infinites- 
imal flexes that there are three surfaces, z = 1, w = 1 and n = 1, in 
^(C^s). Their union in fact forms the entire RUM spectrum of the grid 
framework, since the form of the matrix function (for the single vertex 
motifs) is diagonal with entries 1 — 2, 1 — w, 1 — m. However the point 
we underline here is that a line-localised periodic flex leads directly to 
a hyperplane of wave vectors in the RUM spectrum. 

Similar observations hold in three dimensions for plane-localised flexes. 
In three dimensions, for example, such a flex, which is assumed to be 
"in-plane phase-periodic", leads to a line of RUM wave vectors. This 
in the case for the framework for perovskite in Example (h) below. 

In the case of exotic RUMs a diversity of examples are implied by 
recent experimental and simulation studies in material science. (See 
in particular Dove et al We now note that the matrix function 
perspective also leads to the expectation that the set of RUMs of a 
"typical" material crystal is likely to be exotic. 

In view of Theorem 15.11 Proposition 13.11 and Deflnitions 12. 5[ 14.11 the 
low energy vibrational RUM spectrum of a crystal framework formed 
by a 3d zeolite will be approximated by the zero set, on the 3-torus, 
of the three- variable polynomial pc{^) of the framework. If the edge 
motif set F^. has n edges with an "external" vertex not in then this 
polynomial is "likely" to have degree n. At one extreme, this fails if 
there exist localised flexes, since then the polynomial vanishes identi- 
cally and at another extreme, in the absence of symmetry for example, 
the polynomial may fail to have any zeros on the torus, apart from the 
point (1, 1, 1) which is ever-present in view of the acoustic/translational 
modes. However, away from these extremes there will be generally be 
RUMs arising from the surfaces and curves of intersection deflned by 
irreducible factors of Pciz). Linear factors create planes and lines and 
(roughly speaking) one expects that for complex crystals that there will 
not be a full linear factorisation of the crystal polynomial and that some 
of the nonlinear factors will contribute RUMs. This was conflrmed by 
Wegner [22] in the case of the framework for the HP tridymite phase 
of Si02 for example. We give another explicit exotic factorisation below 
for a two-dimensional zeolite. 
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5.3. The RUM spectra of crystal frameworks. The foUowing def- 
inition is convenient. 

Definition 5.2. A mathematical zeolite in two (resp. three) dimen- 
sions is a crystal framework C in the plane (resp. M^) consisting of 
congruent triangles (resp. congruent tetrahedra), each pair of which 
intersect disjointly or at a common vertex, or edge, or face, and such 
that every vertex is shared by two triangles (resp. tetrahedra). 

One can view the frameworks of material and mathematical zeo- 
lites at the databases [23] , [23] respectively. Of particular interest are 
the material zeolites, which are denoted by an upper case triple, such 
as FAU, PAU, SOD, RWY, ... , these often being a mnemonic for a 
material group indicator (faujasite, paulingite, sodalite, . . . , etc.). We 
write CxYZ ioi the mathematical zeolite framework associated with the 
material zeolite XYZ when the geometry is unequivocal, as is the case 
for the cubic form of SOD for example (see below). Note that while 
the databases often highlight the framework formed by "T atoms" at 
the tetrahedral centres (each being 4-coordinated with neighbour T 
atoms), we identify the tetrahedral net framework, corresponding to O 
atoms as vertices. 

Example (a): The regular 4-ring framework. A fundamental 2d 
zeolite Ca is determined by the four-pointed star indicated in Figure 
4. It is an instructive example in that it is sufficiently simple that 
one can deduce its RUM spectrum and also its polynomial PCai^i''^) 
from infinitesimal arguments, by the band-limited arguments outlined 
in Section 5.2. 

Note that there are four non-internal edges in the motif and so four 
rows of ^ci^j'w) 'will carry monomials. Thus pc{z,w) has total degree 
at most 4. One can show that there are infinitesimal fiexes of the finite 
sub frameworks as indicated in Figures 3 and 4. The top and bottom 
vertices of each are fixed by the fiex. Each flex extends by horizontal 
periodicity to a band-limited infinitesimal fiex. In the former case there 
is two step horizontal periodicity while in the latter case there is strict 
horizontal periodicity although the band is two cells wide. 




Figure 3. A 2-cen periodic band-limited fiex of Ca- 
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From the discussion in Section 5.2 the first band-hmited flex shows that 
(-l,a;2) e Q{C) for all U2 G T. By symmetry {ui, -1) G Q{C) for all 
Ui G T. Thus Q{C) contains the set 

({1} X T) U (T X {!}) U ({-1} X T) U (T X {-!}). 

It follows that pc{z, w) is divisible by the irreducible factors z — l,w — 
1, z + l,w + 1. Since pc{z, w) has total degree at most 4 it follows that 
either p vanished identically or 

Pciz,w) = (z- l){w - l){z + l){w + l). 

Finally, note that former case does not hold. One can see this, themati- 
cally, by demonstrating that there are no local flexes. Alternatively one 
may simply compute det$(l/3, 1/3) 7^ 0. Thus the RUM spectrum is 
precisely the fourfold union above. 

Examples (b), (c): Deformation phases of C^. The regular 4- 
ring framework Ca has a flnite flex deformation, or phase transition, to 
a companion frameworks Cf, in which each 8-ring of triangles borders a 
regular octahedron. We call Cf, the regular 8-ring framework. A motif 
for the natural translation group is formed by an 8-ring of 24 edges with 
four boundary vertices omitted. See Figure 5. There are 8 edges with 
external vertices and so the expectation is that pc^{z,w) has degree 8 
and this is the case. In particular it follows, from Corollary 14.31 for 
example, that there are no localised RUMs anywhere in the inflnite 
framework. No apparent band-limited inflnitesimal flexes are present 
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in this framework and indeed numerical calculation shows that there 
are no lines of wave vectors and that the RUM spectrum is exotic, 
consisting of four closed disjoint curves on meeting at (1, 1, 1). 




Figure 5. Deformations of Ca, with angle 107r/12 for Cc, 
angle 97r/12 for Cb (in an offset position), angle 87r/12 for 
Ca, rotated by 7r/4, and angle 77r/12, for Cb, the regular 
8-ring framework. 



While the matrix-function ^Cbi^) is rather large, with 24 rows and 
columns, it is sparse and its entries may be conveniently normalised 
by dividing each row by the magnitude of the a;-coordinate difference. 
The magnitudes of the nonzero nonunit entries of the resulting matrix 
are then the tangents of the angles kir /24, for k = 1,3, 5, 7, 9, 11, all of 
which lie in the field extension Q(\/2, y/S). It can be shown that 

PCt{z, w) = z^w'^ + z'^w^ - 2V3{z^w''^ + z'^w + w^z + zw) + 

(4^3 - A){z^w'^ + z'^w^ - 2z'^w^ + z^w + zw'^) + w"^ + z^ . 
and that this polynomial admits an explicit factorisation as a product 

PCb(^, W) = Pl{z, W)p2{z, W) 

where 

Pi{z, w) = z^w — {\/?) + \/2)zw'^ + z + w, 

P2{z, w) = Z^W — (V3 - V2)zw'^ + 2{V3- V2- l)zw -{V3-V2) z + w. 
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Each of the factors is responsible for two of the four closed curves 
that comprise f2(Cb) which is illustrated in Figure 6. 




Figure 6. The curved wave vector spectrum of Cb. 



This approach leads to a rigourous algebraic determination of the 
exotic RUM spectrum but it is nevertheless of interest to obtain a more 
intuitive understanding of such curves and their modes. One approach 
to this is in terms of RUM spectrum evolution which we now sketch. 

Figure 5 indicates four frameworks, say Ci,C2,C3,C4, determined 
by their unit cells. These frameworks are snapshots in a continuous 
smooth path, or deformation, t ^ Ct from Ci = Cc to C4 = Cb which is 
driven by the indicated angle. The angle starts at 107r/12 and decreases 
to 77r/12, the other snapshots being for 97r/12 and 87r/12. Notice that 
C2 is a diagonal translate of and that C3 is a rotation of Ca by 7r/4. 
These frameworks are endowed with the natural isometry group gener- 
ated by horizontal and vertical motions, with differing period vectors. 

We have seen the RUM spectrum of C3 in Example (a). However, 
here the unit cell is doubled and is rotated by 7r/4. Accordingly the 
new RUM spectrum is determined as the image of the old one under 
the "doubling map" {wi, W2) — )■ (w^i, followed by axial rotation. It 
follows that 

VL{C^) = {{w, w), (w, -w):we T}. 

In terms of reduced wave vectors this corresponds to the subset of the 
unit square [0, 1)^ given as the union of the two diagonals. 

On the other hand one can show, by infinitesimal analysis, or by 
numerical calculation that ^(Ci) is standard and, in terms of wave 
vectors, is simply the subset of the unit square [0, 1)^ given as the 
union of the axes. We can now view the RUM spectrum of the regular 
8-ring framework C2 as one picture in a path of sets from the axial 
set fl{Ci) to the cross diagonals set ^(Cs). Indeed the path of spectra 
t — > ^l{Ct), which we refer to as a RUM evolution, is the path of zero 
sets of a path of polynomials and can be analysed in such terms. 
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Example (d): A 2d zeolite with order N. Figure 7 shows the unit 
cell for a 2d zeolite in which there are six framework vertices on the 
boundary. One can verify that there is a nonzero infinitesimal flex of 
the enclosed finite framework with zero displacements at the boundary 
vertices and a vertical displacement at the central vertex. Thus the 
framework has a local RUM and the RUM spectrum is all of T^. 

In general it need not be the case that an order N crystal framework 
has a local RUM internal to the unit cell. This is the case for the 
Roman tile framework for example. Also here one could take a new 
unit cell in which the central vertex is shifted to the boundary and in 
this case one has to consider a threefold supercell before a local RUM 
appears. 




Figure 7. Unit cell for a 2d zeolite with order N. 



Numerical calculation shows that the mode multiplicity function takes 
the following values on T x T; 

iy{l, 1) = 4, z/(l, w) = 3, for 7^ 1, z/(w;, 1) = 2, for w ^ 1, 

iy{—l,w) = 2, and h'{wi,W2) = 1 otherwise. 

In particular the RUM spectrum is standard and so also amenable to 
analysis in terms of band-limited flexes. 

Example (e): Roman tiling framework. The next figure shows 
part of a two-dimensional framework Croih which derives from a classi- 
cal Roman tiling of the plane by regular triangles, squares and hexagons. 
Here the triangles are the rigid units of the framework. Also the frame- 
work is 4-regular and so in particular CRom is in Maxwell counting 
equilibrium. The matrix function for a motif associated with a unit 
cell parallelogram enclosing six vertices is a 12 by 12 matrix in two 
variables. 

The framework evidently has many infinitesimal flexes. Indeed for 
each linear square-hexagon channel there is an (edge-length-preserving) 
smooth deformation {pK,k{t)) in which the channel undergoes parallel 
collapse. Such a motion gives an infinitesimal flex u = [p'^ ^(0)). More 
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Figure 8. Roman tiling framework, with order A^. 

apt to our discussion is the strictly periodic (phase (1, 1)) infinitesimal 
fiex in which each triangle undergoes identical infinitesimal rotation 
about its centre. 

Less evident perhaps are the localised infinitesimal fiexes which are 
supported in a finite region. Indeed there is an interesting minimal such 
fiex u which is zero outside a 2-fold supercell (of four parallelograms) 
and which imparts an infinitesimal clockwise "hexagonal motion" on a 
ring of six triangles in this supercell. It follows that Crohi is of order N 
with full RUM spectrum and that the crystal polynomial is identically 
zero. 

Curiously, despite the multi-deformability of the framework, the 
mode multiplicity function fi : fi(CRom) — ^ is simply given by 
/i(l, 1) = 2 and /i(cj) = 1 for a; 7^ (1)1)- Equivalently, in terms of 
the rigid unit mode wave excitations indicated in the next section, for 
the zero wave vector there is simply the two-dimensional space of triv- 
ial translation (acoustic) modes, while for each other wave vector there 
is a single rigid unit mode. This mode, say for wave vector k = {ki, /C2) 
in [0, 1) X [0, 1), corresponds to the infinitesimal fiex v derived from u 
in terms of the well-defined ("slightly overlapping") infinite sum 

(P,<7)ez2 

Example (f): /3-cristobalite. The so-called kagome net framework 
in three dimensions may be constructed in a layered manner as follows. 
Place upward tetrahedra on alternate triangles of the two-dimensional 
kagome framework Ckag and place downward tetrahedra on the other 
triangles to create a framework layer. Such layers are then fused to- 
gether to fill space, creating the crystal framework we denote as Cknet- 
Note first that the band-limited periodic infinitesimal fiex arguments 
above can be applied to the two-dimensional kagome framework, just 
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as in Example (a), to arrive at its polynomial directly; 

Pkag{z,w) = {Z- 1){W - l){z - w) . 

Note that a phase-periodic flex of the 2d framework, with phase (ci;i, c<;2) 
say, induces " layer-limited" infinitesimal flexes of the kagome net frame- 
work which are phase-periodic. It follows from the discussion in 
Section 5.2 that for all G T there is a rigid unit mode u of Cknct with 
phase {ui,U2,u}). Similar assertions hold for the other axial planes. 
Any natural motif giving rise to the crystal polynomial pknet{z,w,u) 
shows that it has total degree at most 6. Since our arguments show 
that it must vanish on the six planes z — 1 = 0,w — 1 = 0,u — 1 = 
0,z — w = 0,w — u = 0,z — u = Oit follows that 

Pknet{z, W, U) = {Z - l){w - l){u - l){z - w){w - u){z - u) , 

for the monomial order with z > w > u. 

The polynomial above was also obtained in [22] and [16] by cal- 
culation. The framework is also the tetrahedral net framework for 
/9-cristobalite, [1] (and also happens to be associated with the recent 
material HerbertSmithite). 

We remark that the kagome net has other deformations, or phases, 
and one such phase features in tridymite which was the first system 
where curved surfaces of RUMs were observed [1] . 

Note that any 2d zeolite C which admits a two-colouring, with no 
like coloured triangles adjacent, similarly determines a 3d zeolite V by 
the layer construction. In particular this applies to Ca and its phases 
CbyCc- Once again since the layers are effectively decoupled a RUM 
phase {u2, 0U2) for C leads to a line {u2, 0J2, oj) in the RUM spectrum of 
V. 

Example (g): Sodalite. We now define the sodalite framework Csod 
which is a zeolite framework in three dimensions, with cubic symme- 
try, which is built from 4-rings of tetrahedra. (The tetrahedra repre- 
sent Si04 rigid units with silicon atoms at the centres and oxygen atoms 
at the vertices.) These 4-rings, so-called secondary building units in 
zeolite assembly [23], are oriented in the high symmetry arrangement 
indicated in Figure 9. Six such rings may be placed on the the six faces 
of an imaginary cube so that vertices sit on the midpoints of the edges 
of the cube. This gives a finite bar-joint framework that we call the 
sodalite cage. With unit edge length for the tetrahedra the cube has 
sidelength 1 -|- \/2 while the unit cell for Csod has side-length 2 + \/2. 

A motif for the framework can be based on the set Fg of edges in 
three pairwise-connected 4-rings of the sodalite cage. The images of the 
edges of Fe under the action of the associated isometry group T, with 
orthogonal period vectors of length 2 + a/2 are disjoint and generate 
the crystal framework Csod- For the set F^ one could take the vertices 
Fe except for the 12 vertices furthest from the imaginary cube faces. 
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There are 24 such edges and the resuhing matrix function has 72 rows 
and 72 columns. 




Figure 9. The top 4-ring of the sodahte cage. 

We prove that Csod is of order A^. Specifically we show, by infin- 
itesimal flex geometry, that there is a nonzero infinitesimal flex v of 
the finite sodalite cage framework such that all the outer vertices are 
fixed by v. That is, Wk,<5 = if ^^,5 is an outer vertex of the cage. The 
"outer fixed" sodalite cage framework has 36 free vertices with 108 de- 
grees of freedom while there are 144 constraining edges. Despite this 
considerable over-constraint there is sufficient symmetry to allow in a 
one-dimensional space of infinitesimal flexes. 

We show below that an individual 4-ring, Ri say, of the sodalite 
cage has an infinitesimal flex, v^^"^ say, which fixes the quadruple of 
outer vertices (the upper vertices in Figure 9.) In this flex the mirror 
quadruple vertices have equal magnitude displacement vectors which 
are directed towards two opposing corners of the imaginary cube. See 
Figure 9. Taking v^^^ so that these vectors have magnitude 1 it follows 
that v^^"^ is determined up to sign and that this sign may be specified by 
labeling the cube corners "a" and "r" for their attracting and repelling 
sense. Note that one can label the eight corners of the imaginary cube 
in this manner so that no like labels are adjacent. In this case the 
implied flexes v^^\ . . . , w*^^-* of the six 4 rings of the sodalite cage have 
equal displacement vectors at common vertices. This consistency show 
that there is an infinitesimal flex of the entire sodalite cage in which 
the outer vertices are fixed, as required. 

It remains only to show that there is the stated flex of the 4-ring Ri 
One can confirm this by computation or, more revealingly and themat- 
ically, by the following simple symmetric infinitesimal flex argument. 

Let pi,p2 be two non-opposite top vertices of Ri with intermedi- 
ate vertex ps, let pi,P3,P4 be the vertices of an inward facing face of 
a tetrahedron of Ri with vertices Pi,P3,P4,P5, so that the lower ver- 
tex p5 is a cube-edge midpoint. There is a unique (oblique) "upward" 
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displacement velocity of the intermediate vertex which has unit 
length and is such {ui,U2,Us) = (0,0,^3) is a flex for the two edges 
[Pi7P3]7 [^37^2]- The displacement vector M3 induces a unique displace- 
ment vector M5 which is in the direction of the cube edge and is such 
that 



The triple ui = 0,us and M5 now determine the inflnitesimal motion 
of the tetrahedron, with flex vector U4 for p^. However, the reversal 
(sign change) of M3 induces the reversal of M5 so it is clear from the 
symmetric position of the tetrahedron that M4 must be the unique unit 
norm (oblique) downward flex at p^. Continuing around the ring it 
follows that Ri has the desired inflnitesimal flex. 

One can apply similar arguments to other zeolite frameworks and of 
course to any such frameworks that contain a sodalite cage. Also we 
note, as do Kapko et al [12], that Crwy is derived from Csod simply 
by replacing each tetrahedron by a rigid unit of four tetrahedra. Thus 
the same inflnitesimal flex geometry applies and Crwy has order N. 

Example (h): Perovskite. Let Cqct be the crystal framework 
formed by corner linked octahedra in maximal symmetry arrangement. 
More precisely, Cqct has the integer translation group T and motif 
{Fy,Fe) where 



= {0, 1/2, 1/2), (1/2, 0, 1/2), (1/2, 1/2, 0)} = K,o : 1 < «: < 3} 



and Fe consists of the twelve framework edges between the centres of 
adjacent faces of the unit cell [0, 1)"^. Thus the vertices of V{Fe)\Fy 
have the form pi^^^, where 71 = (1,0,0), 72 = (0,1,0), 73 = (0,0,1). 

The framework is not in Maxwell counting equilibrium, being over- 
constrained, but has considerable symmetry. The frameworks feature 
in cubic perovskites, such as SiT03, and RUM distributions have been 
determined experimentally [5], [1]. To derive Q{Coct) analytically by 
flex analysis we flrst obtain Q{Csq), since the square lattice appears in 
two dimensional cross-sections for the three axial directions. 

Performing row operations on $Csq7 as given in Section 2, we see that 
{'z,w) is a point of the RUM spectrum if and only if the equivalent 
matrix 



(uS - U3,P5 -P3,) = 0. 
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has rank less than 4. This occurs if and only if 
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has rank equal to or 1. The rank is if and only ii z = w = 1, corre- 
sponding to the two-dimensional space of rigid motions (translational 
infinitesimal flexes) with phase (1,1), and the rank is 1 if and only if 
z = w = —1. Thus 

0(Csq) = {(1,1), (-1,-1)}. 

The infinitesimal flex for the phase (—1,-1) is the one for which the 
rigid units, in this case squares with diagonals, rotate infinitesimally in 
alternating senses. 

Note now that this alternating rotation flex of Cgq induces a plane- 
localised flex of CocT in each of the framework planes x = 1/2,?/ = 
1/2, z = 1/2. By the argument in Section 5 it follows that f2(CocT) 
contains the three sets of phases, 

Tx{-I}x{-1}, {-l}xTx{-l}. {-1} X {-1} X T. 

That the spectrum is no more than the union of these sets and the 
singleton (1, 1, 1) can be seen from a row analysis of the 12 by 9 function 
matrix $Coct(-^) same style as the argument for Cgq. Thus, in 

wave vector formalism, the RUM spectrum of the octahedral net Cqct 
is the set of lines 

(«,l/2,l/2), (1/2, a, 1/2), (1/2, 1/2, a) 

together with the acoustic wave vector (0, 0, 0). 
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